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Abstract 



We perform a model-independent analysis of the data on branching ratios and CP asymmetries of B — > 
<j)K and B — ► ?/ modes. The present data is encouraging to look for indirect evidences of physics 



beyond the Standard Model. We investigate the parameter spaces for different possible Lorentz structures 
of the new physics four- Fermi interaction. It is shown that if one takes the data at la confidence level, only 
one particular Lorentz structure is allowed. Possible consequences for the B s system are also discussed. 
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I Introduction 

O 

The data from B factories has reached a state where one can talk about looking for indirect signals of physics 
beyond the Standard Model (BSM) [1, 2]. Before the Large Hadron Collider (LHC) comes up, it seems we need 
,-SZh ' a fare share of luck to directly produce BSM particles, be it of supersymmetry, technicolour, compactified extra 
dimensions, models with extra chiral/vector fermions, gauge bosons or scalars, or of some other extensions of 
the Standard Model (SM). On the other hand, radiative corrections induced by such new particles may affect 
D . low-energy observables. 

The effects of specific models of BSM on low-energy data, particularly data from B factories, have been exten- 
sively discussed in the literature. Indeed, there is some justification for such an enthusiasm. It seems that the 
data is not exactly in harmony with the SM expectations. There are two caveats: the error bars are still large 
enough, and the low-energy dynamics is not very well known. The second point can be tackled by looking at 
the mixing-induced CP asymmetry data, where hadronic uncertainties mostly cancel out. Unfortunately, only 
for a limited number of channels such asymmetries have been measured. It is hoped that in near future we may 
have a situation where the error in the data will be dominated by theoretical uncertainties. We can only hope 
that the experimental numbers that we now have will stand the test of time and proceed accordingly. 

If one looks carefully, the anomalous results seem to be concentrated in the b — > sss decay sector. One can 
say that this is apparently so because we do not really understand the penguin dynamics, but the radiative 
decay b — > S7 is well explained within the SM. This has prompted a number of analyses of the nonleptonic 
data, particularly those of B — > 4>K, in the context of specific new physics models, like different versions of 
supersymmetry [3], and non-supersymmetric models like extra Z 1 ', more Higgs doublets, extra fermions, etc. 
[4]. A good guide about the anomaly is the Heavy Flavor Averaging Group (HFAG) website [1] (also the UTfit 
collaboration [2] performed some model- independent new physics analyses). The value of sin(2/3), as extracted 
from all charmonium modes, is 0.725 ± 0.037 (this is consistent with the value obtained from CKM fits without 
superimposing the direct CP asymmetry measurement data), while it drops down to 0.43 ± 0.07 for b — ► sss 
modes. Both sets arc internally consistent as far as exclusive modes are concerned. However, there is a second 
important factor: the branching ratios (BR) for B — ► rfK and B — » -qK* channels arc abnormally large, at 
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least if we do not take into account the a posteriori explanations of a large charm content of r\ or rf or the 
significant gluonic contributions. Some of the papers in [3] also deal with this BR enhancement problem. The 
data is shown in table 1. 

Still, the results are not inconsistent with the SM. If we look at table 1, it is seen that sin(2/3) extracted from 
B — > <j)Ks is consistent with that extracted from the charmonium modes at less than 2a. The BR data is never 
a clear indicator; e.g., in the factorisation model [5] the B <j>K BR is quite unstable with the variation of N c , 
the effective number of colours, and one can always take help from SM dynamics beyond naive valence quark 
model to explain the BRs of r](') modes. Only the value of sin(2/3) from B — > r/'K is away from the charmonium 
value by almost 3<t. It is imperative to have more precise values of the CP asymmetries shown in table 1, and 
also accurate CP asymmetries from channels like B — > 3Ks or B — > K + K~Ks [6]. 

Let us assume that the data, as it is, will stand the test of time. Even then this docs not mean that only 
the b — ► sss channel is affected by BSM physics. Indeed, it will be very difficult to find a model which affects 
only one single channel. In particular, if at least two of the final state quarks are left-handed, then the SU(2) 
conjugate channel, b — ► ccs, will be there. The deciding factor is the SM contribution: for the former, it is a 
penguin contribution, whereas for the latter, it is a tree-level one, so the BSM amplitude has less chance to 
compete and show up in b — > ccs decays. 

In this letter we perform a model-independent study of all b — > sss channels that have shown signs of anomaly. 
There are a number of works in the literature [7, 8] that stress the need for such a model-independent analysis. 
We investigate what type of Lorentz structure of the four- Fermi interaction coming from the BSM physics can 
simultaneously solve all the anomalies. For simplicity, we do not include tensor currents, but all other bilinear 
covariants are taken into account. This serves as a testing ground for BSM models as far as their ability in 
explaining the anomalous B decay results are concerned. We will also assume that the new physics, whatever 
it is, does not contribute significantly to B° — B° mixing; in particular, to the CP asymmetry prediction from 
the B° — B° box diagram (e.g., minimal supersymmetry with perfect alignment, i.e., all flavour-changing 5 
parameters set to zero). 

We will show that it is impossible to find a simultaneous solution for all BRs and CP asymmetries with a 
hierarchical structure of new physics (i.e., only one Lorentz structure is numerically significant). The reason 
has been explained in section 3. We find that it is necessary to have some dynamics within SM but beyond the 
naive valence quark model (NVQM) to explain the BRs of the 77C) modes. Even with such an assumption, we 
have a fairly constrained parameter space, different for structurally different four-Fermi interactions, which can 
be further tested in hadronic machines (now that BTeV has been cancelled, this essentially means LHC-B). 

Theoretically, a new Lorentz structure as discussed above should contribute to B s — B s mixing. This can 
potentially change the weak phase coming from the B s — B s box from the SM prediction of almost zero. 
This will directly affect the measured CP asymmeries in B s — > J/tpcf>, B s — ► DfD~, B s — > rf ^rf ' channels. 
Unfortunately, the allowed parameter space for new physics can contribute only marginally in B s — B s mixing; 
the SM contribution, whose lower limit only exists, is too large. Thus, the effects can only be studied (hopefully) 
in a supcr-B factory, an e+e - collider sitting on the T(55) resonance. 

The first two decay channels, i.e., B s — > J/ip4> and B s — > D+D~ will be further affected from BSM physics in 
decay if an SU(2) partner b — > ccs channel is present. However, we will show that this effect can be neglected. 
On the other hand, the channels B s — > <jxj> or B s — ► rf ^rf \ controlled by b — > sss transition, should show 
modification from the SM prediction on two grounds: BSM physics in mixing as well as in decay. These 
predictions will be quantitatively different for different Lorentz structures. 



II Data, Theory, and New Physics 

The relevant data, after Moriond 2005, is taken from the HFAG website [1] (updated April 2005) and is shown 
in table 1. 
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Tabic 1: Experimental data on b — + sss modes. For our convention of Acp, see text. The third column shows 
the value of sin(2/3) as extracted from mixing-induced CP asymmetry of the corresponding modes. The error 
bars are at la confidence limit (CL); the upper limits are at 90% CL. 



The theoretical numbers, particularly those for the BRs, have some inherent uncertainties that stem from low- 
energy QCD. We use the numbers from the conventional factorisation model with N c = 3 [5, 9]. However, one 
may question the validity of such a simplified approach. Other models like QCD factorisation [10] or perturbative 
QCD [11, 12] differ in their predictions of the BRs for these modes. The reason for this difference is mainly 
twofold: first, the form factors for the factorizable diagrams are different, and second, the nonfactorizable 
topologies (like annihilation or emission) arc given unequal importance in these models. A study shows that the 
predictions for the amplitudes vary at most by 20% (i.e., more than 40% variation in the BRs) for the iV c -stable 
modes like r](')K, while 4>K modes are rather unstable and the amplitude uncertainty can go up to 40%. Since 
there is no concensus about which model one should use, we take the conventional factorisation numbers with 
an error margin of 20% in the amplitude level (and 40% for (f)K modes). This is just to take into account the 
variation of predictions for different models and should not be confused with the variation of a model parameter, 
say, A^ c . The charm contribution in r\ or 77' is neglected to start with, since that is an a posteriori approach to 
fix the BRs with the experiment. However, we will soon see that one must entertain either the possibility of a 
significant charm content in these mesons or a large gluonic contribution, i.e., something beyond the NVQM, 
to satisfy the data even in the presence of new physics. 

The CP asymmetries, particularly the mixing-induced ones that measure sin(2/?), are more or less free from 
theoretical uncertainties. We assume that charm and up penguins are negligible compared to the top penguin, so 
that in the SM, B — > 4>K$ is essentially a one-amplitude process. (Again, one can assume some phenomcnological 
values for these penguins which reproduce the CP asymmetry data, thereby alleviating the need for any new 
physics.) The CP asymmetries are defined as 

= T(B+/B^f)-T(B-/W^f) 
CP T(B+/B° -» /) + T(B-/B° -» /) ' 1 ' 

which differs in sign from the HFAG convention. More precisely, defining 

" (f\B d ) ' (2) 

where <f>M is the mixing phase ((f) m = 2/3 for B° — B° mixing in the SM), we have 

d 1-|A| 2 2ImX 

and the conventional S and C(A) parameters are given by 

S = -a% P , C = -A = a d CP . (4) 
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For numerical analysis, we follow the values given in [5]; e.g., the decay constants (in GeV) are 



f n = 0.133, f K = 0.158, U = 0.233, (5) 

and we use f\ = 1.10/ ff , fs = 1-34/^ with the one-angle mixing scheme (the two-angle one [13] makes little 
difference) 

/ 8 cos 6* _ /i sin g ./gcos^ _ /i sin 6> 

_ / 8 sin0 /iCOs6» o /8sin0 /i cos 9 

Here 6* = —22°, and if%p% = (ri\uj^(l — 7 5 )w|0), and similarly for other decay constants. The numerical values 
are = 0.099, /* = -0.103, ,Q = 0.051, /*, = 0.133. Note that /* and are both large and of the same 
order of magnitude; this will be relevant for our future analysis. 

The expressions for the SM amplitudes can be found in [5]. The regularisation scale is taken to be /i = rrib/2 w 2.5 
GeV and the quark masses according to [5] (m s = 0.5 GeV, m u = nrid = 0.2 GeV). We use the Wilson coefficients 
evaluated for iV c = 3 at next-to-leading order (NLO) . Considering the fact that we have considered some inherent 
uncertainty in the amplitudes, it is expected that the results will be more or less stable with the variation of 
N c . This is indeed found to be the case, but for more comments the reader is referred to the next section. 

We discuss two different types of effective four-Fermi interactions coming from new physics: 

1. Scalar : C new = h s e llps (s a (ci + c 2 j 5 )s a ) (sp(c 3 + 0475)6/3) , 

2. Vector : £ new = h v e %4 ' v (s a ^{ci + c 2 ^)s a ) (s/37m( c 3 + c^)bp) . (7) 

Here a and (3 are colour indices. The couplings h SiV are effective couplings, of dimension [M]~ 2 , that one 
obtains by integrating out the new physics fields. They are assumed to be real and positive and the weak phase 
information is dumped in the quantities <f> 8iV , which can vary in the range 0-27T. Note that they are effective 
couplings at the weak scale, which one may obtain by incorporating all RG effects to the high-scale values of 
them. The couplings C\-C\ can take any values between —1 and 1; to keep the discussion simple, we will discuss 
only four limiting cases: 



(i) (S + P) x (S + P) [{V + A) x (V + A)} 
{if) (S + P)x(S- P) [{V + A) x (V - A)} 
{iii) (S-P)x(S + P) [{V — A) x (V + A)} 
(iv) (S-P)x(S- P) [{V — A) x (V — A)} 



Ci = 1, c 2 = 1, c 3 = 1, c 4 = 1; 

Ci = 1, c 2 = 1, c 3 = 1, c 4 = -1; 

Ci = 1, c 2 = -1, c 3 = 1, c 4 = 1; 

ci = 1, c 2 = -1, c 3 = 1, c 4 = -1. (8) 



This choice is preferred since the 1 — (+)7s projects out the weak doublet (singlet) quark field. For the doublet 
fields, to maintain gauge invariance, one must have an SU(2) partner interaction, e.g., s(l — 75)5 must be 
accompanied by c(l — 75)0. No such argument holds for the singlet fields. 

We have chosen the interaction in a singlet-singlet form under SU(3) C . The reason is simple: one can always 
make a Fierz transformation to the local operator to get the octet-octet structure. Since no spin-2 mesons are 
involved, the transformations are particularly simple, as we can neglect the tensor terms. 

The amplitudes for the various decay processes under consideration take the following form: 

(i) B — > (f>K (the contributions are same for neutral and charged channels): 



iM s = ih s e L4 >° (-J—^j ^(cic 3 - c 2 c 4 ),4 , 



iM v — ih v e 1 ' 



' ' 2X) ' 2X' jr; 



(9) 
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where 



= 2f <i> m^F Q B ^ K (m 2 )(e <i> . PB ). 



(10) 



(ii) B — ► rf ) K (again, since there is no b — > sww new physics (NP) contribution, charged and neutral channels 
are equally affected): 



iM s = -ih s e^ s 
iM. v = —ih v e 1 ^" 



c 2 c 3 R(P, b, s) - j^- {(c 2 c 3 - cic 4 ) + (c 2 c 3 + cic 4 )i?(F, b, s)} 



C2C3 + (( C 2 C 3 + C1C4) - (c 2 c 3 - cic 4 )i?(P, 6, a)) 



-4, 



where 



R(P, b, s) = mp/2m s (mt — m s ) 
(P = 77,??') is the chiral enhancement factor, and 

Ap = if P {m 2 B - m 2 K )F (q 2 ), q = PB-PK- 



(11) 
(12) 

(13) 



(iii) B — > 77^ 'if*: This follows a path similar to that discussed in (ii); however, only the axial-vector current 
contributes, and the expression for Ap looks slightly different: 



Ap = 2m^f P {e*. PB )A Q {q 2 ). 
The chiral enhancement factor is given by S(P, b, s) = —mp/2m s (mi > + m s ). 



(14) 



III Result 




Figure 1: Allowed parameter space for (S + P) x (S + P) (lower red) and (S — P) x (S — P) (upper green) type 
new physics at 2a confidence level. The NP amplitude is to be multiplied by 10 -8 . 

We are now ready to discuss our results. They are shown in figures 1-3 as the allowed parameter spaces for 
different Lorentz structures. The theoretical inputs have been summarised in the previous section, while the 
experimental inputs (BRs and CP asymmetries) are all listed in table 1. Note that since we are interested only 
in the b — ► sss transition, no other decay modes (like B — > 7T7t or B — > ixK) have been taken into consideration. 
Although theoretically b — > ccs modes (like B — > J/ipKs) can be affected by the NP dynamics, the justification 
for not including them in the analysis (which essentially means neglecting the NP amplitude in the decay) comes 
a posteriori from the allowed amplitude of NP; this is so small compared to the tree- level Cabibbo- allowed SM 
amplitude that one can safely disregard them. 



Apart from the BRs and CP asymmetries, we also use the following results from [1]: 
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Figure 2: Allowed parameter space for (S + P) x (S — P) (lower red) and (S — P) x (S + P) (upper green) type 
new physics at 2a confidence level. Also shown is the la level parameter space (upper blue) for (S— P)x (S + P) 
type. The NP amplitude is to be multiplied by 10~ 8 . 




Figure 3: Allowed parameter space for (V + A) x (V + A) (red) and (V + A) x (V — A) (green) type new physics 
at 2<T confidence level. Note that at the bottom left corner, there is a large overlap between these two parameter 
spaces. The NP amplitude is to be multiplied by 10~ 8 . 
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• sin(2/3) from charmonium modes: 0.725 ± 0.037; 

• sin(2/3) from the fit to the Unitarity Triangle plus e K [2]: 0.729 ± 0.042 

• sin(2/3) from b — > sss transitions: 0.43 ± 0.07; 

• 41° < 7 < 68°; dependence on 7, coming through the SM amplitudes of B — > rf ) modes, is completely 
negligible. 

The magnitude of the CKM elements are taken from [14], with only the unitarity constraint imposed. All 
elements except V u b and Vtd are taken to be at their central values, while the full uncertainties for these two 
elements were incorporated. Again, due to the nature of the transition we are considering, this variation has 
almost no effect. 

With the expressions for BRs and CP asymmetries at hand, we perform a scan on the new physics (NP) 
parameters h s ^ and <j) s ( v y If we assume the NVQM to hold for the SM amplitude of B — > no common 
region is found that satisfy all the data. This leads us to relax the BR constraint on the B — » r](') channels, i.e., 
we assume that there is some SM dynamics beyond NVQM that generates a significant fraction of the amplitudes. 
We also assume that this dynamics does not affect the CP asymmetry predictions, i.e., in the absence of NP, 
the mixing-induced CP asymmetry in B — > rfKs should give sin(2/3) and the direct CP asymmetries in all these 
channels should vanish. The implementation of such a logic is performed by demanding that the contribution 
of the SM amplitude (NVQM part) plus that of the NP should not overshoot the experimental BRs (we find 
that the beyond- NVQM dynamics should generate 30-40% of the decay amplitudes for these modes). 

Only the regions that satisfy all constraints are shown in figures 1-3. For the scalar-type new physics, it can 
be seen from figures 1 and 2 that all Lorentz structures are allowed when we take the experimental data at the 
2cr confidence level (CL), and their ranges are also in the same ballpark, but what differs is the weak phase. 
This behaviour can easily be explained from the expressions of the NP amplitudes. Since we do not consider 
SU (2) conjugate channels, there is no way to discriminate, as far as the magnitudes are concerned, between 
SU(2) singlet and 51/(2) doublet type interactions. We have not taken into account any strong phase difference 
between the NP amplitude and that of the SM; since the direct CP asymmetries are compatible with zero for 
almost all channels (except B — > rjK by about la) this is not such a drastic assumption. We have checked that 
if we switch the strong phase on, even the B — > rjK CP asymmetry data can easily be accomodated. However, 
strong phase differences cannot alleviate the necessity of introducing beyond- NVQM dynamics. If we take the 
experimental errors at la instead of 2a, only the (5 — P) x (5 + P) Lorentz structure remains allowed. 

It is easy to understand why a simple NP model, without any contribution from beyond-NVQM dynamics, 
cannot explain the BR anomalies of the rj and rj' modes. If we take only a b — > sss structure, the NP amplitude 
contributes equally to both B — * rf and B — > rj, the only difference being the factors /* and /*,, which happen 
to be numerically close. Thus, if the NP amplitude is taken to be so high as to jack up the B — > rfK BR to the 
experimental value, it also raises the B — > rjK BR, which, however, is expected to be much smaller. The weak 
phase being fixed, we cannot even use the argument of constructive interference with the SM in one case and 
destructive interference in the other. Thus, one is forced to admit the possibility of some SM dynamics beyond 
NVQM to explain these decays. 

For the vector channels, only two Lorentz structures are allowed: (V + A) x (V + A) and (V + A) x (V — A), 
both at 2a. The structure (V — A) x (V + A) is a near miss; one can have a tiny overlap if the data (particularly 
the CP asymmetry data) shifts a bit. Thus, it is not entirely prudent to rule out such a structure right now. 
However, if the experimental error bars are taken at la, all the allowed structures get ruled out. So, if the 
numbers remain as they are, and the error bars get smaller, there will be only one possible Lorentz structure of 
new physics. However, it is too early for such a prediction. 

There are many well-motivated models that can generate such Lorentz structures. For example, minimal Super- 
symmetry can produce such vector-axial vector interactions, so does the model with extra gauge bosons in the 
left-right symmetric models. R-parity violating supersymmetry can produce tree-level scalar type interactions. 
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We have also checked the robustness of our results with the variation of the regularisation scale from \i = nib/ 2 
to fi = nib, and found it to be minimal. This is not unexpected since at the NLO level the scale dependence of 
the Wilson coefficients becomes small. Finally, we have found that with the variation of N c , some regions that 
are disallowed for N c — 3 become allowed, but the result is more or less stable. For the scalar case, the results 
are absolutely iV c -stable (with a slight variation of the parameter space) , when we vary iV c over the range 2 to 
4. All structures are allowed at 2a CL but only one, (S — P) x (S + P), at la CL. For N c = oo, there are no 
allowed regions for any Lorentz structure, but that is only to be expected, since the B — > <j)K NP amplitude 
vanishes at that limit. 

For the vector case, with N c = 2, the (V + A) x (V + A) and the (V + A) x (V— A) parameter spaces that we got 
for N c = 3 remain almost unchanged. The other two structures now get allowed over a tiny parameter space, 
namely, h v < 0.4 x 10~ 8 and 7r < <j> v < 2ir. The reason is that with the lowering of N c , the C2-dependent term in 
eq. (9) becomes more important. These results are true for 2a CL. Only the first two remain at la CL, over the 
range 0.3 x < h v < 1.1 x 10" 8 , 20° <</>„< 100° (for (V + A) x {V + A)), and 1.2 x 10" 8 < h v < 1.6 xlO" 8 , 
28° < <f> v < 40° (for {V + A) x {V — A)). For N c = 4 the results are almost identical with N c = 3. The reason 
for this is essentially the highly unstable nature of the B — > <j)K BR with the variation of N c (remember that 
the variation of N c is not the same as keeping N c — 3 and taking different models of mesonic decays). 

The result is stable with the variation of the other parameters (quark masses, form factors, CKM elements, 
etc.). 

We have not discussed here the polarisation data of B — > VV {e.g., B — > <pK*) channels. This will be treated 
in detail in a subsequent publication [15]. We would just like to point out that a scalar- type new physics may 
significantly reduce the longitudinally polarised fraction in the decay B — > <j)K* [16] (for the experimental data 
see [1]). 



IV Signatures in B s System 

We expect a large number of B s decay channels to be probed in the hadronic machines. The first ones will 
include the gold-plated channel B s — > J/ipcj) an d B s — > D+D~, both occurring through the tree- level decay 
b — ► ccs. In the SM, the mixing- induced CP asymmetry is expected to be very small. It has been pointed out 
that new physics may introduce a nonzero CP asymmetry. There are two sources for that: (i) a new amplitude 
in the B s — B s mixing, with a nonzero weak phase, and (ii) a new contribution in the b — > ccs decay. It may so 
happen that both of them are present, and one needs a careful disentanglement of the two effects [17]. 

It can be shown that the new physics amplitude, given by h s / v \, is at most 10% compared to the corresponding 
SM quantity for b — > ccs decay, namely, (Gf/ V2)V cs V c b- This, as far as the decay is concerned, NP effects will 
be extremely hard, if not outright impossible, to detect. This limits our discussion to NP effects in B s — B s 
mixing only. Note that this is not true if we consider final states like <jxj> or r){')n{') that proceed through b — > sss 
transition. 

Unlike the decay, the calculation of the mixing amplitude depends on the precise structure of the NP Lagrangian, 
and one must perform a model-dependent study. Let us take the la result, i.e., the only possible Lorentz 
structure is {S — P) x {S + P). One of the prime candidates to generate such a structure is R-parity violating 
(RPV) supersymmetry, which leads to a slepton- mediated amplitude proportional to A- 22 A- 23 *. It can easily be 
seen that 

V V * 

h s- • (15) 

Thus, to be consistent with the decay data, the magnitude of the product of these two A's should be somewhere 
between 0.6 x 10~ 3 and 1.2 x 10~ 3 , and the associated weak phase between and it (the minus sign takes 
care of a phase factor of exp(i7r)) for a sneutrino mass of 100 GeV. The necessary formalism to evaluate RPV 
contributions to the neutral meson mixing amplitude has been worked out in [18, 19]. Following these and 
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taking AM S > 14.4 ps 1 one can calculate the effective CP-violation in B s — B s mixing. 

The result, however, is not really encouraging. The reason is simple: the SM amplitude is too large, 
enhanced from the B° — B° mixing amplitude by at least a factor of |Vj s | 2 /|V^| 2 times £ 2 , where £ = 
Ib s \fB~B~J ' fs d \fB~B~ d = 1.21 ± 0.06 [20]. Furthermore, only a lower bound on the SM amplitude exists. The NP 
amplitude with such a weak coupling as obtained from the decay fit cannot compete with the SM amplitude. We 
find sin(2/3 s ), the effective phase from the B s — B s box, to be never greater than 0.1. Thus, it will be extremely 
unlikely to have a signal of this type of NP in B s — > J/ip<j) decay. Note that for such a Lorentz structure, there 
will be a small but nonzero NP contribution to the b — > ccs decay. 

One can, of course, integrate out the left-handed squark field and get the decay B s — > tf£^ (this is a strongly 
model-dependent statement and is true only in the context of RPV supersymmetry) . For squarks at about 
300 GeV, leptonic decay widths of B s will be greatly enhanced; in particular B s — > [i + will be close to the 
experimental limit. (This set of couplings also enhances B — > K(*^£ + £~ decay widths, but the constraints are 
more model-dependent.) Thus, these channels will be worth watching out. 
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